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Abstract 

We present a non mean-field model which undergoes a magne- 
tostriction phase transition in the temperature. That is, the crys- 
tal becomes sharply contracted and magnetized once the temperature 
passes below the critical value. 
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1 Model and main theorem 



Magnetostriction is know in physics as a phenomenon of a drastic change of 
geometric shape of crystals, which is accompanied by magnetic transition, 
see e.g. jK], jM]. Usually it is a first order phase transition with a jump 
of spontaneous magnetization together with the jump in geometry of the 
crystal elementary cells. Physical origin of this phenomenon is related to 
so-called magnetoelastic coupling, i.e. to the interaction between spin and 
displacement degrees of freedom in magnetic crystals, jKj. Various mean- 
field theories of this phenomenon were discussed in literature since a long 
time. See, e.g., iZPj and references therein, for crystals, and GZ for magne- 
tosriction in ferrofluids. (The solvable model with a short-range interaction, 
discussed in jM], does not exhibit the jump specific for magnetostriction, 
because it is one-dimensional.) 

In the present paper we propose a simple - and a first non mean-field! - 
model of this phenomenon. We prove that our model undergoes the phase 
transition, when the crystal becomes sharply contracted and magnetized, 
once the temperature passes below the critical value, provided the dimension 
is at least two. 

We consider the following model: at each site s of Z d we have an Ising 
spin er s , while at each bound I = (st) of the lattice we have positive real 
variable r st , playing the role of the spatial distance between two sites. 

Initially we were interested in the Hamiltonian 

H A (a A , r Ab ) =- J2 J M + V Yl ^ - R f- h Y. a " 

(st)eA b (st)eA b seA 

Here the function J(-) > describes the dependence of the strength of the 
interaction between the spins a s and at on their spatial separation. The 
parameter R is the ground-state distance between sites in the absence of the 
spin interaction; h is the external magnetic field. We were assuming that J is 
small on large distances and large on small distances. Our hope was to show 
that in the symmetric case - h = - the model would undergo the striction 
transition as the temperature goes down. But we were unable to show that, 
and, moreover, our computations suggest that such first order transition does 
not take place for the Hamiltonian H. 

To realize our program we have to modify our Hamiltonian, adding an- 
other "geometric" term to the interaction. Namely, we will consider the 
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model, defined by the following Hamiltonian: 

H A (a\ r A >) =- ]T J{r st ) a s a t + p ]T (r st - i?) 2 (1) 

(st),(s't)£A b :\s-s>\=V2 seA 

Here in addition to the parameter p > 0, which is enforcing the lattice 
structure with the spacing to be close to R, we add another parameter A > 0, 
which has the effect of making the r-lattice more regular. In particular, this 
term makes the "triangle inequality violation" energetically unfavourable. By 
the "triangle inequality violation" we mean, for example, the situation when 
among the four bonds r st ,r s 't,r s ' t >,r st ', forming a plaquette of the lattice, 
there are three relatively small values and one relatively big. 

To ensure that the above model undergoes the striction transition we have 
to suppose that the interaction J is weak enough on large distances r, and 
is strong enough on small distances. Otherwise this function can be fairly 
general. We will describe now one specific choice of the class of interactions 
J, for which the transition takes place; other choices are also possible. 

We are supposing that above some value p > the interaction is weak: 

J (r) < u for r > p, 
with u small. We further suppose that the interaction is bounded: 

max J (r) = U < oo, 

r>0 

and that within the region r < p it is sufficiently strong: for some K C [0, p] 
and for all r G K 

J (r) > U, 

with U large, while jj = 1 + x with x small and mes {K} > p/2. As we show 
below (see (|14p). the choice of the parameters R,p,U,u and x is possible, 
which guarantees the striction transition to happen. 

The Hamiltonian (^J) has the Reflection Positivity (RP) property with 
respect to reflections in the shifted coordinate planes: 

Li-k = {x G W 1 : Xi = k} , i = 1, d, 
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with integer k; it is also RP with respect to reflections in the diagonal planes 



again for k integer. To simplify the computations we will use the latter; this, 
however, is applicable only in 2D case. The general case can also be treated, 
using the RP in coordinate planes, along the same lines. 

To formulate our results, we introduce the indicators of some events: 

for a bond I = st we define 



where p > is a parameter to be chosen later. Similarly, we define the 
indicator 



where e > is another parameter to be chosen. 

We call a Gibbs state (•)* , corresponding to the Hamiltonian (JTJ) and 
inverse temperature (3, a contracted state, iff for every I 



Likewise, we call a Gibbs state an expanded state, iff for every / 



Theorem 1 Let h = 0. It is possible to choose the parameters of the Hamil- 
tonian Q in such a way, that the following holds. 

• at all temperatures low enough there exists a contracted Gibbs state; 

• at all temperatures high enough there exists an expanded Gibbs state; 

• for some critical temperature (3 C there exist at least two different Gibbs 
states, {-) C p c and (-)^ ; the state (-)^™ is contracted, while the state (-)^ 
is expanded; 



L 



'i,j;k 



[x G M d : Xi -xj = k] , i, j = l,...,d, i ^ j, 





(pn e > I 
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• if there exists a contracted state at the temperature then in fact 
there are at least two such states, (-)p and (-)^ . They are oppositely 
magnetized: for every s, t 

3 

Our result makes the following conjectures very plausible: 

• above the critical temperature T c every Gibbs state of our Hamiltonian 
is expanded, having zero magnetization, 

• below T c every Gibbs state is contracted, while every pure state has 
non-zero magnetization, 

• at T = T c precisely three pure states coexist: one is expanded, with 
zero magnetization, while the other two are contracted and oppositely 
magnetized. 



2 Basic estimates and proof of the main re- 
sult 

Our strategy of the proof is to follow the RP theory of the first-order phase 
transitions. To this end we introduce the following indicators: 

Pf - of the event {ri =st < p}, 

P® - of the event {p < ri =st < p + e}, 

Pf - of the event {ri =st > p + e}, 

so + Pj° + Pf" = 1. We also introduce the indicators P^,P^ and P^, 
which are products of the above over all bonds, i.e. P£ = YlieA b e ^ c - ^ e 

put p° A > = n eA6 + pn ■ 

The strategy consists in showing that for the finite volume states (-)^ with 
periodic boundary conditions at inverse temperature (3, uniformly in volume: 

• the expectation (Pf)^ is small at low temperatures, 

• the expectation {Pf)» is small at high temperatures, 
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• the expectation (P^Pj?)^ is small at all temperatures and for all pairs 
of bonds I 7^ I', 

• the expectation (Ff) „ is small at all temperatures. 

The rest then is standard, see jSj. 

(1) First we show that the expectation (Pf*) * is small at low temperatures. 

,p>\ ^ V2|A| 

(p l > h<(PDT lA] < ) Ab 



We have : 

(^^-TgyE / II dn/ {n >^}(n)e-^^ 



< V TT / d ne -^(n-«) 2 +^ < 1 2 |A| J—) e 

z M)y£iU z m W/w 



2|A|/3u 
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Here and in the following we use the identity: e Ax dx = \J~\- On the 
other hand 

< P a% = zlB)J2 [ II dra {n < p} (rOe-^ A ^ (2) 
1 /p «ra- Kj R 2 -2A« 2 )\ 2|A| 



Z A (/3) V2 

Therefore, 



(pn P ^ - — > ( 3 ) 



2 v /2^ e -/3(c / -«-^ 2 -2Ap 2 ) 

which is small for (3 large once 

U -u> /ii? 2 + 2Ap 2 . (4) 

(2) Next we show that the expectation (P{^)p is small at high tempera- 
tures: 



/P<\ < /p<\ 1 /2|A| < I ( P A)p 



1/2|A| 



We have: 



■<li(0) 

For the lower bound we have 



ie\ b 



A r A,) (5) 



2 |A| TT / drie -P(n+8\) ri 2 -f3u 



^a(/5) VV /?(M + 8A) 

where we use in the third line the inequality (x — y) 2 < 2x 2 + 2y 2 , and also 
the fact that for every / the sum Yli'-i nn v ( r i ~ r i') 2 nas 4 terms. Therefore 



< /2/j(/i + 8A)^ ( ^ +M . Mp _ fi)2); 



(6) 

7T 

which is small for small f3. 

(3) Now we estimate the correlation function (P l < Pf > ) /3 , with / = (st), 

I' = ( s 't), \s - s'\ = y/2. We have: 



>\1/|A| 



where the indicator P£ corresponds to the following event: on half of the 
bonds - - of A 6 - namely, on those which have one endpoint on the 
sublattice, generated by the vectors (1,1) and (2, —2) - the event r. > p + e 
happens, while on the remaining ones - = A h \ Af - the event r. < p 
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happens. Therefore 

^) p =zmY, J e"^ A(CTA ' rA6) II II ^/ {n >, +£} (n)^i {n ,< p} (nO 



< 1 olM J3(u+u~p( P ~R) 2 -\e 2 )\A\ |A| f /°° . -^{r-Rf 

~ Z k {(3) 6 ' ' 



|A| 



|A| 

To estimate the partition function from below we note that P£ (■) + 
P° A > (•) < 1, so 

Z A (P) > Z K {ff) ((P^ + <P?> > p ) . (7) 
Using (j2J), ©, we thus have 

|A| 



2 I jn_ oe p(y+ u -^p- R ) ~ Xe 



p £ ) „ < ; —w- ( 8 ) 



By suppressing one of the terms in the denominator of (JHJ) we get the 
following two estimates: 

/p<p>\ <r 8 l~^~ e l3(U+u-2U-p{p-R) 2 +2p,R 2 -Xe 2 +4Xp 2 ) / Q \ 

which is good for /3 large, and 

(i^i?^ < 2 \j^P (V + 8A ) e^ U+3u -^ R)2 - Xe2 \ (10) 

which is good for (3 small. So we have to look for some intermediate value of 
(3*, such that for (3 > (3* the r.h.s. of Q is small, while for (3 < (3* the r.h.s. 
of (fTUj) is small. Of course, such value of the inverse temperature should be 
the one which makes the two terms in the denominator of (JHJ) equal; in other 
words, the reasonable choice of the value (3* is to take it to be the solution 
of the equation 



P J3(U-pR 2 -2Xp 2 ) n ~/3u /-|-,\ 

2 6 "V/^(/i + 8A) e • Uij 
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But any choice of (3* would be as good as this one, provided only that the 
estimates Q and (fTU|) will turn into bounds strong enough. 

(4) The last estimate we need is that for the expectation (Pf) „ . We have 



'Jf>, < (pi) 



0\1/2|A| 



Now 



i r r ,v~^ A(ffA ' rA6) 



ieA b 

p+e 



< 1 o |A| ( / dre -P[n(f>+£-R) 2 -u)y\A\ 
Z\{p) \ J 



>a(/3) 

Combining with the estimate (J7J) we find: 



W) B < ^- ; 1 rw ^ 



V2ee 

2|A| 



(P e /3(U-vR 2 -2Xp 2 )Y m + / ZE e -/?« 



0(M+8A)' 

Here we can proceed as in the previous case, turning (|T2*|) into two different 
estimates, depending on the value of (3. However, the case of the observable 
Pf is easier, and it is sufficient to keep just one summand in the denominator 
of (|12j) in order to get a reasonable estimate on it. Namely, we keep the second 
one, arriving to 



(P?) p < ^ 2 ^^ 8X K e-^^- R ^ 2u ]. (13) 

We now shall show that if we make for the Hamiltonian the following 
choice of the interaction parameters : 

A = fj, = 1, U = 2R 2 , U = (2 + 5 2 ) R 2 , u = 5,p = R-\s = 25R, (14) 

with R big enough and 5 small enough, then the conditions of our theorem 
hold for the interval [/?+,/?_] , provided (3 + = (3 + (R,S) is small enough, and 
/3_ = /3- (R, 6) is large enough. 
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Since the estimate (jlj) is satisfied under our choice (JHJ), the relation (J3J) 
holds for all (3 large enough. As we said before, the r.h.s. of (JBJ) is small for 
all j3 small enough. Therefore it is enough to check that the r.h.s. of (JHJ) and 
(112)1 are small uniformly in all (3. 

To proceed with the estimate of the correlation function (P l < Pf > ) „ , as 
indicated above, we have to choose a value of the intermediate inverse tem- 
perature P*. Our choice is 

(3* = 2R~ 2 \nR. (15) 

One can check that thus defined /3* is indeed an approximate solution to 
(111)1 as R —>■ oo, though this is not important. 

In the region f3 > (3* we will use the estimate Q, which under the choice 
(JT3J) becomes 

(P l < Pf'} < 8R ll^ e i 3 *(( 2 + s2 ) R2 + s ~ 4R2 -( R - R ~ 1 ) 2 + 2R2 -' is2R2 +' iR ~ 2 ) 

V $ 

V P* 



<8 # UH—JT'i 1 **) 
\j 2\nR 

< R- 4S2 

for R large. 

In the region (3 < P* we shall use the estimate (fTUjl . which similarly 

becomes 

< lSy^R-h^ 1 - 252 )* 2 



< 18V27T- 1 \nRR~ AS 

< R~ 3S2 



for R large. 
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Finally we consider the bound (|13|) . which becomes 

< P 0^ < Q^l 6Re ~^R-^28R-Rf-28] 




7T 

Note that the function ^/xe~ ax has its maximum at x — which equals to 
j-. Applying this to the last expression with x = PR 2 , we get 



which is small for small 5 at any (3. 

3 Magnetization 

Here we will prove the last statement of our theorem: the occurrence of 
spontaneous magnetization in the contracted states. To do this we split the 
event {r i=st < p} into four events, and we introduce the corresponding four 
indicators 

P ; <±d= - of the event {ri =st < p,a s = ±, a t = ±}. 

We will show now that for all (3 the expectations (P ; <+ ~)^ = (Pj <_+ )^ 
are small, uniformly in the volume. Together with the obvious statements 
that 

(P, <++ ) e = {!?-), 

and 

<^ <++ >, + (rr-) P + {pr-) + <^- + >, = (p?) P , 

that implies our claim, due to the first part of our theorem and by subsequent 
application of the Theorem XX of [SJ. 
We have 

where h the indicator of the event that for every bond I' we have r;/ < p, 
while 

+ 1 if lsi| + |«2 1 is even, 



-1 otherwise. 
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We denote this spin arrangement by o±. So 



l£A b 

- Z A (f3) \ pe J 
As in (0— QUj) , we have two estimates: 

r . r .—fili(p—R\ 2 . . 

/p<-+\ < P C - n p -P{U+p(p-R) 2 ~f,R 2 -2\p 2 } 



and 

-(3p(p-R) 2 



><-+\ / 



/3( M +8A) e 



In fact, with our choice (JUjl of the parameters the second one is effective for 
all p. We have 

\ * //3 V 7T 

The r.h.s. has its maximum at (3 = 2{(r-r 1 - 1 ) 2 -s] > so 

<P<- + >, < 3iT 
which is small for i? large enough. I 



>-i 



1 



27re [(i? - J?- 1 ) 2 - 5]' 
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